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Abstract 



An integral solution to the quantum Knizhnik-Zamolodchikov (gKZ) 
O ' equation with \q\ = 1 is presented. Upon specialization, it leads to a con- 

. jectural formula for correlation functions of the XXZ model in the gapless 

regime. The validity of this conjecture is verified in special cases, including 
the nearest neighbor correlator with an arbitrary coupling constant, and 
general correlators in the XXX and XY limits. 



1 Introduction 



^ ■ Consider the one-dimensional spin 1/2 XXZ chain 



1 oo 

^=-0 E + + A«+l) • (1-1) 

n=— oo 



In this paper we address the problem of describing correlation functions of (|1 . 1| ) 
in the gapless regime |A| < 1. In the earlier works ||l|, the case of the anti- 
ferromagnetic regime A < — 1 was treated in the framework of representation 
theory of the quantum affine algebra Uq{sl2)- As a result, correlation functions 
have been described by using the quantum Knizhnik-Zamolodchikov (gKZ) equa- 
tion. It is this aspect that we will be concerned with in this paper. Before coming 
to the content, let us first recall some known results for A < —1. 

Let V = C^, and consider the R matrix R{(3) G Endc(V^® V) associated with 
the XXZ model (see (p.2|)). The qKZ equation is the following system of linear 



*Department of Mathematics, Faculty of Science, Kyoto University, Kyoto 606, Japan. 
^Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606, Japan. 



1 



difference equations for an unknown function • ■ ■ , P2n) that takes values 

in \/®2n. 

GniPl, - 2ni, ■ ■ ■ ,/?2n) 

= Rjj+i{Pj - Pj+i - 27ri)~^ ■ ■ ■ Rj2n{Pj - P2n " 27ri)"^ 

xi?i,(/5i - Pj) ■ ■ ■ Rj-ijif3j-i - f3j)Gnif3u ■ ■ ■ , /^i, ■ ■ ■ , P2n)- (1.2) 

Here Rij{(3) G Endc(^®^'^) signifies the matrix acting as R{(3) on the (i,j)- 
th tensor components and as identity elsewhere. The correlation functions of 
arbitrary local operators are obtained as the specialization 

n n 

Gn(P + 7rt,---,P + m,P^"^). (1.3) 

To be precise, in the case A < —1, there are two functions F^*^ (i = 0, 1) 
associated with the two anti-ferromagnetic vacuum states, and it is their sum 
Gn = Fj^^^ + Fji^^ that satisfies the gKZ equation ( |1.2| ), as well as a set of relations 
(|2.4|) , (|2.5| ) and ( p. 61) . The correlators are given by specializations of F^'^ rather 
than Gn itself. In the context of representation theory, the functions F^'-* are 
traces of products of certain intertwiners (vertex operators) taken over the inte- 
grable highest weight modules ^(A,). By realizing V{Ai) in terms of bosonic free 
fields, an explicit integral formula was obtained for the functions F^*^ {i = 0, 1) 
and hence for the solution Gn of the qKZ equation. 

The argument relating correlation functions to the functions Gn is based on 
the extension of the corner transfer matrix method f^, |^ . It is applicable to the 
more general case of the XYZ spin chain as well [^. Correlation functions are 
related in the same way as above with solutions Gn of the qKZ equation, this 
time having the elliptic R matrix as coefficients. Unfortunately the mathematical 
structure of the XYZ model is not fully understood yet (see ||^, ^ for a formulation 
of an elliptic extension of Uq{sl2))- The free field realization is still unavailable 
(see however the recent development J^, in this direction). Thus it remains 
an important open problem to construct solutions to the gKZ equation in the 
elliptic case. 

For the XXZ chain in the gap less regime |A| < 1, the corner transfer matrix 
fails to be well defined. Nevertheless this case can be viewed as a limiting case of 
the XYZ chain, so that the same recipe ([l73|) is expected to apply for obtaining 
correlation functions. (Unlike the case A < — 1, the vacuum state is unique and 
the distinction between F'-"^ and F*^^-* disappears.) The problem is then to find 
appropriate solutions of the gKZ equation. 

Up to an overall scalar, the R matrix R{(3) of the XXZ chain is a rational 
function in 

C = e-''^ q = -e-'\ (1.4) 

where 

A = — cosvrz/ = ^ . (1-5) 
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However the nature of the solutions is quite different depending on whether A < 
— 1 or |A| < 1. In the case A < —1, we have — 1 < g < and the solutions are 
meromorphic in (, typically involving infinite products of the form nj^Li(l~CQ'^")- 
On the other hand, the case | A| < 1 corresponds to |g| = 1. There are no analytic 
solutions which are single valued in (. Instead one has to look for solutions which 
are meromorphic in log^. 

A certain class of solutions to the gKZ equation with |g| = 1 has been studied 
in detail by Smirnov in connection with the form factors in the sine-Gordon 
theory. The equation relevant to the correlation functions of the XXZ model is 
slightly different from Smirnov's, in particular the shift —27ri in ( p..2|) is replaced 
by +27ri in his case (the former has 'level —4' while the latter has 'level 0'; we 
will consider also the case in which the shift — 27ri is replaced by —iX where 
A > is a parameter.) In this paper we give a solution to the former, in the 
form of an integral which has a similar structure as in the case A < — 1, and 
conjecture that its specialization (|1.3|) gives the correlation functions of the XXZ 
model with |A| < 1. Our integral formula is essentially the same as the one 
written down earlier by Lukyanov [jlO[- However, in Lukyanov's case, it is given 
as a generating function of the form factors of local operators in the sine-Gordon 
theory. Our point here is to interpret it as a formula for the correlation functions 
on the lattice. In general, difference equations determine the solutions only up to 
arbitrary periodic functions. We need to ensure that the particular solution we 
present actually corresponds to correlation functions. As supporting evidences, 
we verify this statement in three special cases for which exact results are available: 
(i) the nearest neighbour correlation (cxfcrf), (ii) the XXX model A = — 1 and 
(iii) the XY model A = 0. The integral formula and these verifications are the 
main results of this paper. 

The text is organized as follows. In Section 2 we formulate the gKZ and allied 
equations. We then write down an integral formula for solutions. In Section 3 we 
specialize the formula in Section 2 and propose that it gives correlation functions. 
In the special case z/ = 0, we recover the formula for the correlation functions of 



the XXX model derived earlier in O, 12,0]. In Section 4 we process our integral 



formula to reproduce the simplest correlation function (crj^cTg). This quantity 
can be derived by differentiating the ground state energy of the Hamiltonian. In 
Section 5 we consider the XY limit {u = 1/2), which can be studied independently 
by using free fermions. The correlation functions are given by the determinants 
of certain matrices whose entries are elementary functions in f]j. Our integral 
formula in this case is shown to be equivalent to the free fermion result. The last 



Section 6 is devoted to a discussion concerning some previous works |T^, |12 
on the gKZ equation with |g| = 1. 

Since most of the statements are proved by purely computational means, we 
put together technical points in the appendices in order to ease the reading. In 
Appendix A, a summary of Barnes' multiple gamma functions is offered. Ap- 
pendix B contains the proof of the difference equations of Section 2. In Appendix 
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C it is shown how the n-fold integral is reduced to an (n — l)-fold one by exphcitly 
carrying out the integration once. Appendix D is the derivation of the expression 
for ((Ti(t|). Appendix E is the evaluation of an integral in the case v = 1/2. 
Finally, Appendix F is devoted to the free fermion theory in the XY limit. 

Acknowledgement. We thank Sergei Lukyanov, Atsushi Nakayashiki and Feodor 
Smirnov for valuable discussions. We wish to express our sorrow at the death of 
Claude Itzykson. We have always liked his nice lectures given in his fascinating 
voice ever since we first met him in San Francisco, 1979. 



2 Integral fromula 

2.1 The difference equations 

In this section, we formulate the system of equations we are going to study, 
including the qKZ equation with |g| = 1. We then give a particular solution in 
the form of an integral. Throughout this section we fix parameters u and A (see 
(|1.4| )) such that < z/ < 1 and A > 0. For the convergence of the integral, we 
assume that 

X + ^>2n. (2.1) 

In the application to the XXZ model, we will choose A = 27r. 

Consider the /2-matrix R{l3) G Endc(^® V) acting on the tensor product of 
V = Cv+ ®Cv-: 

£i,e2 

R{P) = -^RiP). (2.2) 
The parameter u enters the matrix elements as follows. 

i?++(/3) = :R::(/5) = i, 

Rt-W) = R~tiP) = m, 
Rll{P) = RZm = m. 

R!;if^{(3) = in the other cases, (2.3) 

where 

j,^, sinhz//3 sinhz/vri 

= • ^ • m ' = ■ u ( ■ ^• 

smhu [in — p) smhu [in — p) 

The function K{j3) will be specified below. It is chosen to ensure that the i?-matrix 
satisfies the unitarity and the crossing symmetry relations 

i?i2(/3)/?2i(-/3) = id, Rf:U/3) = RZHiT^i - /3). 
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Let n be a non-negative integer. Consider a V^^^-valued function G„ = 
GniPi, ■ ■ ■ , P2n), depending on the 'spectral parameters' ■ ■ ■ , P2n- We set 

Go = l. 

We study the following system of difference equations for Gn involving the pa- 
rameter A: 



(2.4) 



GniPl, ■ ■ ■ , l32n-l,P2n " '^•^)ei,-,e2n " G^n(/^2n, A, " ' ' ; P2n-l)e2„,ei,-,e2„-iy (2-5) 



Gn{Pl, ■ ■ ■ , P2n) 



£i,---,£2r. 



/32n=/32n-l+'ri 



Gn-l{Pl, ■ ■ ■ , P2n-2) 



(2.6) 



In particular, the qKZ equation (|1.2| ) is a consequence of (p.4|) and (|2.5|) . It can 
be shown also that (|2.4|) and (|2.6D imply 



G'n(/5l, ■ ■ ■ , /32n)£i,--,, 



£2)7 



/3j+l=/3j+7r« 

(2.7) 

for any j = 1, ■ ■ ■ ,2n — 1. Note that the equations ( p.4|) -( pl6D involve only 
the functions G„(/5i, ■ ■ ■ , P2n)ei,--,E2n with fixed value of the 'spin' ei + ■ ■ ■ + e^n- 
Throughout this paper we will restrict ourselves to the 'spin 0' case, i.e. we 
assume 

G„(/3i, ■ ■ ■ , /32n)ei,-,e2„ = UuleSS Ex A h £2n = 0. 



2.2 Auxiliary functions 



Our aim in this section is to construct a solution to ( p.4| ), ( |2.5| ) and ( |2.6[ ) by using 



an n-fold integral. The formula involves certain special functions k{(3), p{(3), (p{(3)^ip{(3). 
Let us first give their definitions and list some of their properties. In what follows, 
Sr{x\u!i, ■ ■ ■ ,ujr) will dcuotc the multiple sine function (see appendix for the 
definition). 



S,m'2^,^)S2{^f-^P\2n,^) 
S2{-tP\27r,^)S2in + zP\2n,^) 

dt sin ^tsinh(l - v)t 
sinh t cosh ut 




= exp 
k{(3)k{-(3) = 1, 

smh ulm — pj 
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n(B) - sinh ./3 '^3(7^-^/3)>g3(7^ + A + ^/3) (o(^^o (^in^ 
piP) - smh -X Ss{-iP)Ss{X + iP) ^^'^''^ - '^^^'^1^''' 

p(/:^j-p(^jexpy^ i cosh sinh i sinh ^ 

^^^^ -.(/^), 



= ^^^^m + TTi) + • • • when ^ -ni. 
i sin TTi/ 



^2(1 + i/3|A4)52(f-i/3|A,^)' 

f /-^ disin^ ^sinh(l + 

ip{(3 - i\) _ sinhi/(/3 - f ) 



ip{l3) sinhi/(/3 + f -iA)' 

V{P±^) _ sinhf(/?±f) 



<^(/?) sinhf(/3Tf ±^) 

'■iSM\l){(^T^i) 



4 sinh 



= sinh f 52(7r + ^)S2{7i - t(3\X, ^) 

+ iX) sinh + iX — ni) 
i/jiP) ~ sinh + TTi) 
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In addition we will use a constant c„ given by 

u 



In the case A = 27r the formulas simplify to 



n(n—l) 

= sinh/3, c„-^ ' 



2.3 Integral formula 



Let us present the integral formula for Gn{Pi, ■ ' ' y P2n)ei-e2n- Given a set of 
indices £1, • • ■ , ^ {+, — }, we define a map a G {1, ■ ■ ■ , n} — * a G {1, ■ • • , 2n} 
in such a way that = + and a < 6 if a < 6. Define further a meromorphic 
function 



Qn(a|/?)ei...e2n 



Uj<a sinh i/(aa - + f ) ni>a sinh z/(/?j - + f ) 



na<i, sinh i/(aa - - ni) 
where a = (ai, ■ ■ ■ , a„) and /? = ■ ■ ■ , /?2n)- After these preparations, we set 

j<k 



X 



n/ ^n¥'K-/?.)n^K-«^)Q"(«i/?)™- (2.18) 



Clearly we have, for any 7, 

GnlA + 7, ■ ■ ■ , ^2n + 7) = G'„(A, ■ ■ ■ , /?2n). 

In Appendix B, we will prove that with the appropriate choice of the inte- 
gration contours (1 < a < n) as given below, the function Gn{Pi, ■ ■ ■ , P2n) is 
meromorphic and satisfies ( |2.4D , (p.5|) and ( |2.6D . 

In order to specify the integration contours, let us examine the poles of the 
integrand of ( p.l^ ). The poles of ip{aa — Pj) are at 

a, = ±z(niA + iii7r + f) {n,,n2>0). (2.19) 

The poles of — «&) (a. < b) are at 

aa-ab = ±iini\ + Mtt) (m, 722 > 1), (2.20) 

and the poles of ■ , , ^ r are at 

aa-at = ^^m (n G Z). (2.21) 
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Since ipi^^a — c^b) has zeros at 

aa - at = ±z(niA + ^vr + tt) (ni,n2>0). (2.22) 
the poles of . /(°''-"'') are at 

aa - ab = i{niX + ^2^71) (ni,n2>l), (2.23) 

or 

«a - Ob = -«(mA + Mtt) (ni>0,n2>l). (2.24) 
Therefore, the poles in the variable of the integrand are contained in the set 

{Pj ± z(niA + mvr + f ) (1 < j < 2n; ni, > 0); 
ab + i{niX + ^^^tt) (a < 6 < n; ni, n2 > 1); 
ttf, - i(niA + H3_:^7r) (a < 6 < ni > 0, ris > 1); 
Ob + i(niA + ^^tt) (1 < 6 < a;rai > 0,^2 > 1); 
Ob - iiniX + ^tt) (1 < 6 < a;ni,n2 > 1)}. 

We choose the contour Ca for by the following rule: 

tta lies on the real line for \aa\ ^ 0, (2.25) 
Pj + f (1 < J < 2n), at + i{\ + i^vr) (a < 6 < n), 

«fe + i^^TT {I < b < a) are above Ca, (2.26) 
/3j - f (1 < J < 2n), a;, - ^i^vr (a < 6 < n), 

afe - z(A + I^tt) (1 < 6 < a) are below C„. (2.27) 



772 



Note that we can choose Ca to be the same contour C for all a such that Pj + 
(1 < j < 2n) are above C and Pj — ^ (1 < J < 2n) are below C. 

Let check the convergence of the integral. Recall that the periods of the double 
sine function 5*2 used in ip and are such that cui = A > and 0^2 = ^ > vr. In 
the proof below, we use "const. " to mean different constants which appear in 
the estimates. From ( [A. 14 ), we have 



\ip{aa — Pj)\ < const, e 



LJT1J2 I "I 



tp{aa - ab) 



sinh ^{aa — — vri 
I sinh ^{aa — Pj ± ^ 
Collecting these estimates we see that 



< const, e ^1^2 d^-l+l""!), 



sinhi/(aa - A ± ^)| < const, e'^;'"'"'. 



■k(2tx—u>\ — 0^2) 

^[pLa - (yb)Qn{(^\P)ei-e2n < COUSt. 6 ^2 

a,j a<b 



Y[ ^{aa - Pj) n V'(«a - ab)Qnio;\P)^^. 



Since we have assumed that 2TT — uji — uj2 = 2tt — X — ^ < 0, the integral is 
convergent. 
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2.4 One-time integration 



In the case of interest A = 2tt, the ra-fold integral for Gn{Pi, ■ ■ ■ , P2n) can be 
reduced to an {n — l)-fold integral by carrying out the integration once. The 
result is stated as follows. 



l<j<k<2n 

71 i 



X ■ 



x-f, /o^ — ^ X^(-l)'^^ /■ • 7 n dakD{ai, ■ ■ ■ , tti-i, ■■■,"«) 
z=i J J k^i 



2n 

X n [n fi^k - Pj) n ^inh //(a^ - /3j + f ) n slnh + + f) 

k^l i=l j<k j>k 

sinh z/ {j2k^i ttfc + I A" - I Ejyr /^j + 7rz(r - 2Z + i)) 
X — : — - — -. . (2.28) 

Ur<s,r,sj^i smh u{ar -as- -m) 

Here we have set 

Cn = 23"(«-i)/2(-vr2p(7r2))-", 

and 

D{xi, Xn-i) = det (e-("-2^-^)^^ 

As before, the numbers 1 < ■ ■ ■ < n are determined by 

{[ I 1 < / < n} = {j I 1 < J < 2n, Sj = +}. 

The integration is taken along a path going from — oo to +oo in such a way that 
—7t/2 < lm{ak — Pj) < TT /2 for all k,j. In the above, we assume that < u < 1/2 
for the convergence of the integral. It should be possible to treat also the case 
l/2<z/<lby introducing a suitable regularization as in but we do not go 
into this question here. 

The derivation of (p.28|) will be given in Appendix C. 



l<j,fc<n-l 



3 Correlation functions 

We now proceed to the description of correlation functions of the XXZ model. 
From now on, we assume that A = 27r. 

First let us set up the notation. Let E^,,! denote the 2x2 matrix with 1 at 
the {e, e')-th place and elsewhere. Thus the Pauli spin operators read 

In the tensor product ■ ■ ■ ®Vj ® V^+i Cg) ■ ■ ■ of l^- ~ C^, we let a", E^^) denote 
respectively the operators acting as cr° or Ef^^i on the j-th component and as 
identity elsewhere. 
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By a local operator we mean an element of the algebra generated by ct°'s. Any 

local operator is a linear combination of operators of the form O = eI^I E^J',^'^] ■ ■ ■ Eifl 

(r < s). The correlation function of O is its expected value with respect to the 
ground state eigenvector of the XXZ Hamiltonian. We conjecture that it is given 
by the following special value of Gn [n = s — r + 1): 



(e: 



(r) 



E^flf^G^(P + m,---,P + 7ri,P,---,P)- 



We shall consider a slightly more general object 



(3.1) 



(3.2) 



This corresponds to introducing spectral parameters f3j as inhomogeneity of the 
model (in the terminology of [jl3l, the corresponding model is 'Z-invariant'). We 
shall denote ( |3.2D by 

(3.3) 

To see the specialization ( |3.2| ) is well-defined, we note the following. In the general 
formula ( p.l8|) the contour C for integration is such that Pj + i{27ini + + 5) 
(ni,n2 > 0) (resp., Pj —i{27cni + ^n2 + ^) {ni,n2 > 0)) is above (resp., below) C. 
When Pj = Pk + TTz, the contour C is pinched by Pj — ^ and Pk + However, 
Qn{c(\P) has a zero at aa = Pj — ^ for d > j, and also at aa = Pk + ^ for a < k. 
Therefore, if j < k, there is no pinching by poles of the total integrand. 



In order for (|3.1| ) to make sense as a correlator, we must check the following 
property: 



Proposition 3.1 



(41 
((E« 



^{E^:l+E^'i)){^, 

,){Pr,---,Ps 



+ E^:l)Ei:^'l^^...E^;lm,pr^i 



Ps-1, Ps) 
■,Ps) 



Proof. Let us take Pj = Pj^i + m in (|2.4| ). Since 



,Ps). 



(3.4) 
(3.5) 



/O 



R{ni) 



we find from (|2.7| ) that 

Gn~l{Pl, ■ 
= T.Gn{P 



, Pj-1, Pj+2, ■ ■ ■ , P2n) 



1:' ' ' : Pj: Pj+1: 



P2n)ei, 



,£2n 



/3,=/3,-+i+iTi 
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Eq.( |3.4| ) is a direct consequence of this. The above equation together with 
imply (since A = 2tt) that 



,£2n-i,e 



/32n=/3i+7ri 



From this follows (p.5|) . □ 

We now write down the integral formula for (|3.2| ). Set n = s — r + 1, and 
define i,---,n(r<l<---<?7,<s + n = r + 2r;, — l)by the following rule: 

{i, - ■ ■ ,n} = {] \ r < J < s, e'j = -} U {j* \ r < j < s, ej = +} 

where j* = 2s + 1 — j . We have then the following expression for ( p.2| ). 

(4l---41)(/?^r--,/5.) 

_ j-j- sinh(/?j - /• /"-Q^ -Q sinh(a,-a,/) 



r<i<fc<s sinh - A) ^ ^ 27r sinh u{ai - a, - vr?) 

^ n [n - f3, + ^0) n -nh.(az-/?.) n sinh.(-a. + /?,+vr. 

r<l<s J=r ^ ' ^ r<]<l l<J<s 

X n r?-— f— m n sinhi.(-az + /3,) n sinhKaz-/3,+ 



Here the symbol a; — Pj + iO (resp. ai — [3j — iO) indicates that the contour for 
ai runs above (resp. below) jSj. 

Let us put the formula in a form closer to the known result for z/ = 0, taking 
r = 1, s = n. We choose the integration contour C+ for Ua {1 < a < n) and 
for (n + 1 < a < 2n) in such a way that (3j + vri (resp., (1 < j < n) are 
above (resp., below) C+ and [3j (resp., [3j — ttz) (1 < j < ri) are above (resp., 
below) C-. (Here the contours are directed from —oo to cx), as opposed to the 

in page 122-3 of ref.||.) 

Set A = {j I e'j = -} and A = {j \ ej = +}. We suppose that \^{A')+\^{A) = n 
since otherwise {eI}]^ ■ ■ ■ E^?lj{l3u ■ ■ ■ , /3„) = 0. We define a mapping 

a e {!,■■■ ,n} = A+L\ A_ ^ a e {!,■■■ ,n} (3.7) 

by the condition that (i) {a | a G A^} = A', {a | a G A_} = A; (ii) if a, 6 G A^ 
and a < b then d <b; (iii) if a, 6 G A_ and a < 6 then d >b. In other words, the 
+'s in the sequence —e[, ■ ■ ■ , —e'^, e„, ■ ■ ■ , ei are — e'j, ■ ■ ■ , —e'g,, eg, ■ ■ ■ , where 
s = tl(^-) = n — s'. Then we have 
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X 



n 



SA. 



n 



smh.{f3j - f3k) 



dan 



n 



sinh(a;a — a^) 



n 



sinh i^(Q;a — 
sinh(aa - 



l<a,j<7i 

-p-j- sinh z/(/5j — fta + Tri) -p-j- sinh — (3j + vri) 
sinh z/(/5j - Ji_ sinhi/(aa 



(3.8) 



In the hmit z/ = we recover the integral formula for the XXX correlation 
functions (P, 0,0]). 



4 Nearest neighbour correlator 



In this section we examine the simplest cases of the general formula for the 
correlators proposed in the previous section. 

First consider the case Gi(/3i,/52). Taking n = 1 in the formula ( p^.28D , we 
immediately find the following. 



Proposition 4.1 

G'l(A,/?2)- 



G'l(/?l,/52) 



^ 1 p(/?i-/?2)sinhf(/ji-/j2-7rz) 
2z/ p(vri) sinh — /?2 — Tri) 

In particular, by setting Pi = P2 + T^i, we have 

{E^l) = (E^l) = 1. 

Next let us take n = 2 in (|]2|). 

Proposition 4.2 Assuming < u < 1/2, we have 



(4.1) 



G2(/5i,---,/54) 



++- 



p(0) V(7r^)^ 27rz/2 E -=1 e'/^i 



(4.2) 



X 



r ^ 

/ ciae" ip{a - Pj) sinh z/(-a + p^ + f) sinh z/(-a + /34 + f ) 

sinh z/(-a + + ^) sinh i^(a + /^2-/3i-fe-/?4 _ 3|i) 
- sinh u{a - /?i + ^) sinh z/(a + ft-fe-fe-z^d _ 

r/ie integral is taken along a path from —00 to +00 such that — ^ < Im (a — /9j) < 
^ /or a// j . 
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Here we have used the relation p(0)p(7ri) = — which follows from ( |2.12| ), 

(CT)and^2(7r|27r,^) = v/2. 

Upon specialization ■ ■ ■ , P^) = {P + ni, P + ni, /3, this integral can be 
processed further. After a chain of steps detailed in Appendix 0, we obtain the 
following result. 

Proposition 4.3 We have 

{E^y_E^y_) = G'2(/3 + 7rz,/3 + 7rz,/3, 

1 d f sinhfl -v)t\ 1 , , 
'sinTTi^/ — r^^M+-- (4-3) 



tt"^ simriy dv \ Jo sinhtcoshi^t 

We note that, since both sides are holomorphic with respect to z/ for < Re z/ < 1, 
([4.3|) is valid without the restriction < u < ^. 

We now compare the formulas ( |4.1| ), ( [4. 3D with known answers. For this 
purpose let us quote from the results concerning the XXZ model which are 
relevant to the following discussion. 

The XXZ model for a periodic chain of circumference N is given by the Hamil- 
tonian 

^ = -2 E (^J^jVi + ^l^l+i + A^I^I+i) • (4.4) 



This Hamiltonian ( |4.4] ) is associated with the six vertex model with the Boltz- 
mann weights (|^, eq.(8.8,9)) 

. jj-w ^ . j^ + w . 
a = sm — - — , = sm — - — , c = sm/i. (4.5) 

Denoting by T{w) the transfer matrix of the periodic system with columns, 
we have 

-^logT(«;) =-^—fH + ^cosfi) (4.6) 
dw w=-fi Zsmfix 2 J 

where A is related to via 

A = -cos/i. (4.7) 

The gapless regime |A| < 1 corresponds to fi being real. 

For a local operator O, let (vac|0|vac) denote its ground state average (in the 
limit A^ — > oo). Since in the gapless regime the vacuum |vac) is invariant under 
the + — symmetry, one must have 

(vac|a^|vac) = (vac| (^Jj - si^l) |vac) = 0. 
Along with 1 = (vac| + -Si^l) |vac), this means 



(vacl-E^^+lvac) = (vac|-E'i"'"l |vac) = -. (4.^ 
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Our formula ( ^.11 ) is consistent with this. 

In the limit — ^ oo, the free energy per site / is given by (@],eq.(8.8.17)) 



kT 



log a + 



sinh(/i + w)xsinh('7r — fi)x 



dx. 



(4.9) 



It follows from the relation 
XXZ chain is 

eo = --cos/.-2sm/.- 



oo 2x sinh ttx cosh jjx 

that the ground state energy per site Cq of the 



(4.10) 



Differentiating cq with respect to A, we can obtain the nearest neighbour 
correlator for the operators: 



(vac I cr^ 1 vac) 



deo 



2 deo 



Inserting 

(vaclcr^cTglvac) 



dA sin fi dfi 
into ( [4.10| ), we find the following expression for this quantity: 



4 d 



sm fi 



sin /i dfx 

On the other hand, in view of (14. 81) we have 



sinh(7r — fi)x 
sinh TTX cosh fix 



■dx 



(4.11) 



vac|aiV^|vac) = (vac|(l - 2Ei^l)(l - 2E^^L)\vac) = 4(vac|Ei'i^i^i|vac) - 1. 



(2) 



^(1) 77.(2) 



Therefore, the formula (|4.3|) agrees with (^.llj) with the identification fi = ttu. 



5 The XY limit 

In this section we study the integral formula for the correlation functions ( p.6|) at 
a special value of the parameter u = 1/2. This is the case where the XXZ chain 
reduces to the XY chain A = 0. It is well-known that the XY chain is equiv- 
alent to the two-dimensional Ising model. To be more precise, the XXZ model 
with A = corresponds to the critical Ising model. In this case, diagonalizing 
the transfer matrix in terms of free fermions, one can calculate the correlation 
functions directly in the presence of arbitrary spectral parameters. The diago- 
nalization is worked out in Appendix ^ Here we show that the formulas thus 
obtained give the same result as the integral formula ( |3.6| ). 
We shall consider the function 

given by (|3.3|) . In order to simplify the presentation, we shall take all /5j's to be 
real throughout this section. (This is a matter of convenience and not actually a 
restriction. The final formulas are valid as meromorphic functions in /3/s.) 
A special feature about u = ^ is that ( |3.6| ) becomes a determinant. 
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Proposition 5.1 If y = \, we have 
{E^l---E^et)iPr,---,Ps)= n 2^coshi(/5,-/3,)xdet(4,-,)i<,,fc,<„. (5.1) 

r<j<k<s 

Here n = s — r + 1, the Iki are given for r < I < s by 

J 2tt 

I ^ s ^ 

^ n 2 cosh i (a - n 2 sinh |(a - (3j + tO) ' ^^'^^ 
and /or s + l<l<s + n 

hi = 4z* 

wzt/i the bar denoting the complex conjugate. In \5.2j , the integration contour is 
a line above the real axis, as indicated by the symbol +iO. 

Proof. Specializing the formula ( p.6|) to z/ = 1/2 we find 

n 2cosh^^/-../n^n2^sinh^(a,-a.) 

r<j<k<s ^ '' •' l=\ KV ^ 

' 1 A 1 



2 cosh \ {ai — /3j) _j 2 sinh — Pj + iO) 



r<l<s 3=r 2V ' ^J'' j 



' 1 A 1 



X TT 2i '^^+' TT = TT 

s+i<LJ ;i 2 cosh liai - (3,) ^j:^ 2 sinh - (3, - zO) 
Inserting 

n 2i sinh -{ai - ai>) = i-nin-l)/2^-^^{a^+-+a„) fka,,\ 

KV ^ - . - 

we obtain the right hand side of (|5.1|). □ 

In Appendix ^ we evaluate the integral ( p.2| ) explicitly (see ( |E.2| )). 
We now proceed to the calculation of correlation functions of the fermion 
operators 

We shall consider only monomials consisting of an even number of such operators. 
They are local operators in the sense of Section 3. For instance 

i^mri=\E^T} {m = l),. 
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Clearly the function (O) for a monomial O is (see ( |3.3| )) unless it consists of 
the same number of V^'s and ip*^^- 

The following two propositions will be proved in Appendix 0. 

Proposition 5.2 



Proposition 5.3 



(5.4) 

^ {m = l). (5.5) 
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The formulas ( [5.4| ), ( |5.5| ) give the same results for the corresponding quantities 
( [F.16|) (vac|?/'m'^r|vac) obtained directly by diagonalizing the Hamiltonian (see 
Appendix ^). In general, the multiple correlators of the fermions are given by 
applying Wick's theorem. Since (vaclV'm'V^/lvac) = (vaclV'^V'rlvac) = 0, the result 
is given as a determinant in the same way as (p.3|). Any local operator (i.e., a 



finite linear combination of monomials in cr"'s) can also be written as a linear 
combination of monomials of the fermions. Therefore, we can state 

Proposition 5.4 For an arbitrary local operator O, (O) = (vac|0|vac) holds. 



6 Discussions 

Before concluding the paper, let us touch upon previous works on the qKZ equa- 



tion with |g| = 1. In Smirnov introduced and solved a system of difference 
equations for the form factors of local operators in the sine-Gordon theory. His 
equations are the same as ( p.4|) and ( p. 51) in this paper except that 5* = —R is 
used (see page 29 in |T^) and that A = —27i instead of A > (the case relevant 
to the correlation function is A = 27r). There is a significant difference between 
( |2.6|) in this paper and the equation (16) (page 11) in Smirnov's. The latter 
requires that the solution has a simple pole at the point = /?2n-i + T^h while 
the former requires that the solution is regular there. The physical origin of this 
difference is that in Smirnov's case the poles are the annihilation poles of the form 
factors while in our case (12.61) is the normalization of the correlation functions 



(see Proposition p.l|) . 
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There are other mathematical differences between Smirnov's formula and ours. 
The number of integration is n-fold in our formula in contrast to the {n — l)-fold 
integrals in Smirnov's. Since the integration can be carried out once (see ( p. 28 )), 
this difference is rather superficial. The significant difference is that in Smirnov's 
formula the (ra — l)-fold integral reduces to the determinant of an (n — 1) x (n — 1) 
matrix with entries given by integrals with respect to a single variable. This is 
not the case in our formula (except for = ^ — the case of the XY model). This 
lack of determinantal structure is already noted by Nakayashiki in [|12| , where the 
special case i/ = was studied. 

In a recent paper [|T5|, Smirnov has constructed an affluent family of solutions 
that corresponds to a family of local operators in the sine-Gordon theory. In our 
case, the structure of the total space of solutions is absolutely unknown. In this 
connection, let us mention an open problem: to show that our integrals satisfy 
G„(/3i, ■ ■ ■ , /?2ri)ei,-,e2n = Gn{Pi, ' ' ' , /^2n) -ei -ean • ^ direct Verification seems 
difficult, and we suspect that it should follow from the uniqueness of the solution 
satisfying certain analyticity and asymptotics. 

In fact, the form factors and the correlation functions are closely related. As 
was discussed in |16, 0, in the regime A < — 1, the former are represented by the 



type II vertex operators and the latter by the type I vertex operators. The type I 
vertex operators generate a family of solutions to the form factor equations, and 
vice versa. In the sine-Gordon theory, this viewpoint was explored by Lukyanov 
in Lukyanov has introduced the appropriate commutation relations for 

the vertex operators, and has given a bosonization of the sine-Gordon theory 
with a cut-off parameter. Though we have not checked the details, it seems 
likely that the integral formula for G„ in this paper is derivable from Lukyanov's 
bosonization after taking the cut-off parameter to infinity. 

In the approach of this paper, the role of the quantum affine algebra Uq{sl2) is 
unclear. In the A < — 1 regime, the free energy, the excitation spectrum and the 
correlation functions depend on the spectral parameters I3j through Q = e~'^^K 
In the gapless regime, these quantities are single-valued only in [3j and the period 
^ is lost. What is the implication of this fact in the representation theory? This 
is an interesting question to be asked. 



A Multiple gamma functions 



The multiple gamma and sine functions were introduced by Barnes |T^, |T8|, Shintani[|T£ 
and Kurokawa||20||. Here we follow the notation of [^. In what follows we fix an 



r-tuple of complex numbers a; = (cJi, ■ ■ ■ , uOj). For simplicity we shall assume that 
RetUj > 0. We set n-ui = riiuji + ■ ■ ■ + nrUJr {r = {ni, ■ ■ ■ , fir)), \ui\ = uji + ■ ■ ■ + uJr- 
The multiple gamma and associated functions are defined as follows. 
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Multiple Hurwitz zeta function 



Cr{s,x\ui)= ^ (n-uj + x) ^ (A.l) 

ni ,■ ■ ■ ,nr>0 



Multiple gamma function 



Tr{x\uD = exp(C(0, x\u)) ('= ^) (A.2) 



Multiple digamma function 

)^(x\iiA — 

dx 



ll) r {x\u) — -^logV j.{x\ui) (A. 3) 



Multiple sine function 

Sr{x\u) = Tr{x\ui)-^Tr{\ui\ - x\ui)^-^^'' (A.4) 



Multiple Bernoulli polynomials 



^ -^Br,n{x\ui) (A.5) 



When r = 1, they are related with the ordinary gamma and other functions 



via 

X 



27r 

iIji{x\ui) = —(ilj{—)+\ogu;i), 
Sx{x\bJ^ = 2sin(— ). 

We list here the basic properties of these functions. 

Difference equations Set = • • • , o^j+i, • • • , a;^)- 

X + Ci^ilo;) — Cr(S) ^1^) = Cr--i('5, x|a;(i)), (A. 6) 

rr(a; + a;i|ci;) 1 



rr(x|aj) rr_i(a:|a;(i)) ' 

+ ci;j|a;) 1 
5'^(a;|ci;) 'S'r-i(:x|u;(i)) ' 

+ LOi\ji) - -B^,„(a;|a;) = n5^_i,„_i(x|a;(i)). (A. 9) 



(A.7) 
(A.8) 
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Analyticity As a function of s, Cr{s,x\ui) is continued meromorphically on the 
whole complex plane and is holomorphic except for simple poles at s = 
1, ■ ■ ■ , r. We have 

(—1)" 
[n — 1)1 

Til 

Cr(-n,x|a;) = (-1)'^7 — — ^Sr,n+r(a;|cu) (n > 0), 
[n + ry. 

lim(s - n)C{s, x\uj) = (-1)"-- (r > n > 1). 

s^n (^12. — l)!(r — n)! 

rr{x\ui)~^ is an entire function of x. Tr{x\Lti) is meromorphic with poles at 
X = n- u [rii, ■ ■ ■ ,nr < 0). 

Sr{x\ui) is entire in x when r is odd, and is meromorphic when r is even. 
Its zeroes and poles are given by 

r:odd zeroes 3.t x = n- u (ni, ■ ■ ■ , > 1 or rii, ■ ■ ■ , < 0), 
r:even zeroes at x = n- u (ni, ■ ■ ■ ,nr < 0), 
poles dX X = n- ui (ni, ■ ■ ■ ,nr > 1). 

All zeroes and poles are simple if n ■ u/s do not overlap. 
Integral representations If Re x > then 

e-^*(-t)^-i dt 



Cris,x\ui) = — r(l — s) 



cnLi(l-e-'^'*)2vr2' 
1)^- f e-^*log(-t) dt 



i—l)' f 

logr^fxlu;) = 7 1 — Brrix\ui)+ I 

r! Jc 



-l^o/ .^1 A f e-^Hogj-t) dt 



where 7 =Euler's constant and T{x) denotes the ordinary gamma function. 
The contour C is shown in Figure |1|. 

Asymptotic expansion Assume ui, ■ ■ ■ ,uJr > 0. Then for any A > 1 we have 

iogr.(zy = (-ir^E^^^(iog. + 7-Ej)+7C(o,.) 

^1 [n + ry. 

as 2; ^ cxD in the angular domain \Arg{z — x)\ < tt — e, where x > and 
< e < TT. 
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Figure 1: The contour C. 



The case r = 2 is of special interest to us. In this case the following formulas 
hold. 

r sinh (x - ^^^) t dt 
= Jc 2 sinh ^ sinh ^ ^^^^"'^ 2^' (0 < < + 



(A.ll) 



2 

S2{X + iUi\iu) _ 1 

S2{X\UD " i^^' 

7TX 7TX 

S2{x\ui)S2{—x\ui) = —4 sin — sin — , (A. 12) 

UJl UJ2 

S2{x\ui) = ^^x + 0{x^) {x — ^0), (A.13) 



i^l+UJ2. , . 
\UJ = 1. 



S2{uJi\ui) = ^j^, 5*2(^1^^^) = V2, 5*2 ( 
In addition, as x 00 (±Imx > 0), we have 

1 Q f \ \ ^ ■ I 1^1+ UJ2 1 CUi u;2 A , . 

log52(x|^) = ±-Kl X - — (— + — + 3) + 0(1), 

\2uJiUJ2 lijJ\ijJ2 12 U)2 tUi ) 

log52(a + x\^S2{a -x\^ = ^ ~ + o(l). (A. 14) 

B Proof of difference equations for Gn 

We prove here that the integral formula ( |2.18|) possesses the required properties 



Proof of Let = GJp^^ p„ = Uj<kp{Pj - Pk)- Because of (|3), (0) 

reduces to th same equation for Gn wherein R is replaced by R. 
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There are four cases to consider: {ej,ej^i) = (—,—), (+,—), (— , +) and 
(+,+)• 

• Case {€j,€j+i) = (-, -) 
We are to show that 

Gn{- ■ ■ , Pj+i, Pj, ■ ■ ■) = Gn{- ■ ■ , Pj, Pj+i, ■ ■ ■) (B.l) 

This is obvious because the integrand of Gn is symmetric with respect to Pj and 
Pj+i if = (~) 

• Case {ej,ej+i) = (-,+) 

Suppose that a = j+l, and set a = aa- Comparing the integrands of G„(- ■ ■ , Pj+i, Pj, ■ ■ ■)...+_... 
and (?„(■■■, Pj, Pj+i, ■ ■ ■)...±=p..., we see that the desired equahty follows from 

b{pj - pj+i) sinh z/(a - Pj + f ) + c{pj - pj+i) sinh iy{pj+i - a + ^) 
= sinh iy{Pj - a + 

• Case {ej,ej+i) = (+, -) 

This is similar to the case {ej,ej+i) = (— , +). 

• Case {ej,ej+i) = (+,+) 
We are to show that 

Gn{- ■ ■ , Pj+i, Pj, ■ ■ ■)•••++■•• = Gn{- ■ ■ , Pj, Pj+i, ■ ■ ■)•••++•■■• (B-2) 

Suppose that a = j, and set a = and a' = Oa+i- Apart from the factors that 
are symmetric with respect to Pj and and antisymmetric with respect to a 
and a', the integrand of the LHS of (|B.2| ) contains 

sinh iy{a' - pj+i + ^) sinh ujpj - a + ^) 
sinh ula — a' — ni) 

Antisymmetrizing it with respect to the variables a and a', we obtain an expres- 
sion that is symmetric with respect to Pj and Pj-^-i- Therefore, we have ( |B.2| ). 

Proof of ^2.3i ). Because of ( p.9|) , the equality (|2.5| ) is equivalent to 

GniPl, ' ' ' , P2n—1, P2n ^-^)ei---e2n ~ G„i^P2n, Pi, ' ' ' , P2n-l) £2n£\---£2n-i • 

(B.4) 

If = —, then n ^ In. In this case, the integrands of the LHS and the RHS 
coincide because of ( |2.11| ): 

V2(aa-/52n+«A) sinhz/(/52n-aa + ^-«A) = V9(aa-/52n) sinhi/(aa-/?2n + ^)- (B.5) 
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If £2n = +, then n = 2n. We make the following change of integration variables: 

iX in the LHS, 



(B.6) 



f ai - 




a2 - 

< 






an- 



Then the integrands become the same by virtue of 
ip{an — — Pj) sinh z/(a„ — iX — [3j + ^) = ^{oin — Pj) sinh — a„ + ^) (B.7) 



and (g:T4D,(gl6D, 



an 



sinh ultta — ttn + iX — ni) sinh z/(q;„ — — ni) 



(B.8) 



We must also check that the contours for the LHS and the RHS are the same. 
Consider the contour Ca corresponding to aa except for the case when e2n = + 
and a = n. (We use Ca and Ca to distinguish the contours before and after 
the change of variables.) The condition (|2.25|) , (|2.26| ) and (|2.27|) for j ^ 2n are 
unchanged for either £2n = + or e2n = —, and for both LHS and RHS. As for the 
case e2n = — and j = 2n, the condition is that, in the LHS, 



/52„ - iA ± i(niA + ^vr + f ) (ni,n2>0) 
are above (below) Ca', in the RHS, 

/32n ± ^(riiA + ^2^7r + f ) (ni,n2>0) 



(B.9) 



(B.IO) 



are above (below) Ca- They are not the same, but not contradictory, i.e., no points 
are required to be in opposite sides of a contour at the same time. Because we 
know that the integrands are the same, it means that those points that appear 
only in either ( |B.9D or ( [B.IOD , are actually not poles. Therefore, the difference 
between ( [B.9| ) and ( |B.10| ) makes no difference in the integrals. As for the case 
e2n = + and j = 2n, the conditions ( |2.26D and (|2.271 ) are unchanged for 
(a ^ n). 

If £2n = +, the contour for is such that for j ^ 2n, in the LHS 
/3j+iA±i(niA + ^7r + 5) (ni,n2>0) 
are above (below) the contour for a„; in the RHS 

/?j ±i(niA + ^7r + f) (ni,n2>0) 
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are above (below) the contour for a„. These two conditions are not contradictory 
in the same sense as above. For j = 2n, the conditions (|2.26|) and ( p.27|) are 
unchanged. 

If £2n = +, the mutual position of aa and a„ changes from the original one 
because of the change of variables. However, the resulting positions of aa and a„ 
in the LHS and the RHS are identical. Therefore, the intagrals are the same. 

Proof of ^TSj ). The factor p{P2n-i — P2n) has a zero at = P2n-i + vri: we see 
from (^) that 

p(/32„-l - /32„) = -^^{(32n - /32n-l - VTz) + • • ■ . (B.ll) 
Sm TTZ/ 

The integral may have a pole at (32n = P2n-i + because the contour Ca is 
pinched by the pole of (p{aa — P2n~i) at cta = P2n-i + ^ and that of ip{aa — P2n) 
at aa = P2n — ^- We will check if this is indeed a pole, and, if so, compute the 
residue. 

Let us consider the four cases separately. 

• Case {e2n-i,£2n) = (-,-) 

The pole of ip{aa — /?2n-i) at = P2n-i + ^ is cancelled by the zero of 
sinh z/(/52n-i — «a + Therefore, there is no pinching in this case. 

• Case ie2n-l,£2n) = (+,+) 

If a 7^ 2n — 1, 2n, for the same reason, there is no pinching of Ca- Consider the 
integrals J, (i = 1,2,3) corresponding to the following contours 

The integral J3 has no pinching at P2n = P2n-i + ^ri. Let us show that Ji — I2 
and I2 — I3 are regular at /?2n = P2n~i + ^ri. After integration with respect to 
«!, ■ ■ ■ , a„_2, the integral reads as 



dan-i [ dar, 

j = 2n — l,2n 

sinh v{an - /?2n-i + y) sinh i/(/?2n - "n-i + f) 



xV'(a„_i - a„) 



sinh v\an-\ — an — vri) 

Here, A{an-i,an) is holomorphic and symmetric with respect to a;„_i and ar. 
Since '(/'(Z^) = we can antisymmetrize the last factor and obtain 



B{an-1, an, (32n-l, hn) 



sinh z/(a„_i - (32n-i + f) sinh v{(32n - a™ + 



1 [ sinh v{an - (32n-i + f) sinh z/(/?2n - "n-i + y) 



2 [ sinh z/(a„_i — a„ — ttz) 



sinh — On-i — ^rz) J 
The integral Ji — J2 is equal to the integral over the contour 
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Figure 2: The contours for Ji, I2 and I^- 



Taking the residue at = /32n — ^ (the minus sign in front of Res below 
comes from the clockwise orientation of the integration contour) , we get 

X(/7(q;„_i - I32n-l)v{0in-l " P2nW{P2n " /32n-l " f ) 

The integrand has no pole at q;„_i = /?2n — ^ because -(/^(a^-i — /32n + ^) vanishes. 
The integral has no pole at /32n = l32n-i + '^i because S(q;„_i, /92n — (^2n-ii P2n) 
vanishes. Therefore, the integral is regular at = P2n-i + Tri. By a similar 
argement we can show that I2 — I3 is regular at (32n — P2n-i + ^ri. 



24 




Figure 3: The contour for Ii ~ I2. 



Case {e2n-l, £2n) 



Taking into account the zero of the factor p{P2n-i — P2n) and the pole of the 
residue —Res 



ni, both at /?2n = /?2n-i + T^i, we have 



C73 



P(/52n-l - /32n)V5(an " /32n-l) (-ReS^^^^^^_^V'(a„ - (32n)) 



Cn-1 

X sinh z/(a„ - (32n-i + ^) 

l<j<2n-2 

X Jl V5(aa - /32n-l)V5(tta - /52n)^(aa - "n) 



X 



l<a<n-l 

sinh z/(/32n-i - + ^) sinh t/(/32n - Qg + ^) 
sinh i^(aa — — vrz) 



1. 



Using ( PTTUQ , ( PT^O , ( P7[B[ ), and Cq = 1, we obtain ( PTZl ). 

The case (e2n-i,£^2n) = (+, — ) is similar. 



C One-time integration 

In this appendix we show how to reduce the n-fold integral for ■ ■ ■ , (32n) to 

an {n — l)-fold integral ( p.28| ). We shall follow the method suggested earlier to us 



by F. Smirnov. A similar calculation has been published in Nakayashiki's paper 
where the limiting case — > was discussed. Since our working is entirely 
similar to the one in |]T2|, we shall only indicate the necessary steps, omitting 



further details. In the sequel we set (3 = {(3i, ■ ■ ■ , (32n), ^ = (^^i, ■ ■ ■ , ^2n)- We 
restrict to A = 2tt, so that ip{f3) = sinh/? in the general formula ( |2.18| ). 
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Step 1 Using 

r>s 

we rewrite the main part of G{P)e as 



k=l k,j 

where 

_ Uj<-k sinh 1/(0;^ - Pj + f ) Uj>k sinh i/(-Q;fe + /3j + f ) 



nf.<s sinh ^(ar — ag — 7ri) 



Step 2 In the first column of the determinant, substitute e ^^"'^ by the right 
hand side of the identity 



■n+lcy2n-l ( n 



Here 

2n 1 ■ 2ra 1 • 

F+(a) = n sinh - (a - + -) , F_ (a) = (-1)" [J «inh -(a- [5, - -) 

and q(/?) denotes some function of /?j's. Terms containing q(/?) vanish in 
the determinant. 

Step 3 Expand the determinant at the first column to obtain 



^n+l22n— 1 n 



=1 



= y'---y'nc^afcn<^K-/5^0(^+M-^-M)A(a)g(Q;|/3),, 

A;=l fej 

where for brevity we set Diipi) — D{ai, ■ ■ ■ , ai-i, ai+i, ■ ■ • , an)- 
Step 4 Next we carry out the integral over a;. For each I, set 

Uj(<i) sinh z/(ai - (3j + f ) ni(>o sinh //(-a; + /3j + ■ 



!!,-(<;) sinh ^{ar — oii — m) Y{r{>i) sinh v{oii — ar — iri) 
Consider the integrals 

2n 
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c. 



C4 





X + vri (s 










R 










X — vri (s 





-R 



,3 



Figure 4: The contour C±. 



taken along the contours C± = Y^t=i C±,i, respectively, which are shown in 
Figure ^ (To fix the idea we draw the figure assuming the f]j are all real, 
but the necessary modification in the general situation should be obvious.) 

It can be verified that inside the contours the only poles of the integrand 
Hi^^ are = =F for s < / or s > / respectively. Collecting the residues 
we obtain 

'(+) _L - _9^„Y A/r(+) 



s{<l) s(>0 



s /' 



Mi? 



Res, 



One can show that, upon integration by the other variables and summing 
over /, these terms cancel with each other. More precisely, set 



<5L(a| 



nj«r) sinh v{ai - (3j + f ) Ilior) sinh u(^-ai + (3j + f 
Y[r«i) sinh z/(ar — ai — vrz) Y[r{>i) sinh z/(q;/ — — ni) 

r^s r^s 



Then we have, for any pair r < s. 



a\ 



+ (-1; 
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This can be seen by changing the variable a,. — + vrz. 
Step 5 From the transformation properties of if{(3) it follows that 



Hi? 



which implies that the integrals corresponding to C-t 3 give the same results 
as for Czp^i. 

As i? — 00, the integrals along C±,2,4 are calculated from the following 
asymptotics of the integrand as ai ±00: 

'^{ai - f3j) ~ 2exp (^^^^^^{ai - f3j)j , 

2n / ^ 2n \ 

QiW)e ~ exp(±z/(naz + AO) X (±1)" 



with 



From the last two steps we find that 
r(+) ui-h,^ - 



where 

Ai = Ai + 

2 . , 

and TZ signifies a term which vanishes when R ^ 00. Hence we arrive at the 



result ( p.28| ) stated in the beginning. 



D Derivation of the nearest neighbour corre- 
lator 

Here we derive the formula ( ^.3| ). We start from ( |4.2| ) and consider the special- 
ization 

G = G{(3 + vri, /5 + vri, /3, /?)++— • 
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Proposition D.l 



2txu'^ J-oo cosh a (p{a) 



(1 — cos TTz/ cosh 2z/a) (ia. (D.l) 



Proof. First let {(3u = {(3 + iri, (3 + iri, f3 + e, (3 + e). Then (|]|) becomes 

—1 e~^"^ piiri — eY 



X 



X 



da e"ip{a — (3 — 'Kif'Lp{a — (3 — eY sinh^ v{a — (3 — e — —) 
sinh u{a — /? —) sinh z/(a — /5 — e —) 



71% 71% 

+ sinh u{a — j3 — —) sinh v{a — (3 — e — —] 

Since a = f3+e+7Ti/2 is not a pole, the contour can be taken as 7r/2 < Im {a— (3) < 
37r/2. Changing the variable a^a + P + e + 7ii and using (f{a + 7!-i)(f{a) = 
—i I cosh a sinh z/(a + 7rz/2) , we find 

G = ,i,„_L^_r<i„e"»'<" + ^'' -1 

e^O 47rZ/2 1 - J-oo 



X 



cosh^a 

7r2 7r^ 7r2 
sinh v{a H he) sinh v[a H ) + sinh via V e) sinh via — 



Noting 







da ■ 



cosh a 



• 1 2 / 7T% -12/ " 

smh + — ) + smh ^[a — — 



711. 



and letting e —>■ 0, we obtain 

1 d / 



+ e)2 -1 



27rz/2 de \J-oo V^(«)^ cosh 



a 



7C% 7C% \ 

xRe sinh u^a + — + £) sinh u^a + ) 
= /1 + /2, 

where /i is given by the right hand side of ( p.l|) and 



£=0 



I, = 



1 r°° . e" _ . . _ , 7T% 



da ., 
47rz/ J-oo cosh a 



2 — Re sinh 2z/(a + — ). 



Using 



'OO 



da- 



7lX 



■00 cosh^ a . ^-^ 
sm — 
2 



(D.2) 
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we find ^ 

thereby completing the proof of the lemma. □ 
Proposition D.2 

G = Ji + J2-i (D.3) 

1 , sinhfl + z/))f: / 1 1 sinhz/t\ 
= tdt H -. Im — + ' 



TT^ 7o sinht ^sinvrz/ cosh z/(t + vri) cosh^z/t^ 



, , sinhfl + z/)t 1 cosnu 

J2 = — : / dt : ^ . I Re 



TTsinTTZ/Jo sinht \^ cosh z/(t + vrz) coshz/ty 

Proof. Substituting the integral formula for log into Ji above, we obtain 

1 _ 1 sinh(l + u)t 

G + — — „ / dt 



2 TT^z/ JO sinh t sinh 2z/t 



X 



/ da -2 — sin ( a ] (cosh 2z/a; cos nu — 1) . 

J-00 cosh a V TT / 



cosh a 

By integrating over a using ( p.2| ), the right hand side becomes 

1 , sinhfl + v)t( / t - Tri t + ni n 2t 

dt — COSTTZ/ ; ; — + * 



vr^ Jo sinhtsinh2z/t y ^cosh z/(t — Tri) cosh z/(t + vri)^ coshz/ty 

After a little algebra we obtain (p.3|) . □ 

Proposition D.3 We have 

^ 1 T T . 1 d f . , sinh(l-z/)t\ 

G-- = Ji + J2-l = T- sinvrz/ / rft — . (D.4) 



2 vr^ sin TTZ/ c?z/ I Jo sinhtcoshz/t 

Proof. Consider the integral 

r sinh(l + z/)(t — 7ii) t — ni ^ 
Jc sinh(t — Tii) cosh ut ' 

where the contour C is as shown in Figure |^. 
Taking the imaginary part, we obtain 

^ , sinhfl + v)t ^ ( 1 \ 
tdt ^ ^Im ' ' 



R sinht ycosh z/(t + vri) y 



+ / dtRe 

JO 



R JCe Je y \ sinh(t — TTz) coshz/t^ 
sinh(l + z/)(i? + ti — vri) R + ti — ni 



sinh(i? + ti — iri) cosh z/(i? + ti) 
sinh(l + v){—R + ti — ni) —R + ti — ni 
smh.{— R + ti — ni) cosh u {— R + ti) 
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-R + ni 



-R 



C 



—e e 



R + TTi 



R 



Figure 5: The contour C and the semi-circle C^. 



As e — > 0, the integral over the semi-circle gives tt^ simny, and as i? — > oo, the 
last term behaves like 



Since 



and 



we have 



Im 



47ri?cos7ri/ — 27r^ siuTTi/ -I- 0(it! 

/ sinh(l + z/)(t — ni) t — 7ii\ 
\ sinh(i — 7ri) coshut J 

t cosh(l + p)t sill Tiv TT siiih(l + u)t cos ttz/ 



sinh t cosh ut sinh t cosh ut 

sinh(l + z/)t sinh z/t cosh(l + z/)t cosht 



sinh t cosh vt 



sinh t cosh vt 



r 

J-R 



sinh t cosh ut ' 
1 



^ , sinh(l + z/)t / sinhi/i 1 

tdt ^ h - — .... , , 

smh t \ cosh i/i sm ttu cosh + ttz) 



-Im 



-f 

J-R 



tdt' 



cosh t 



sinh i cosh^ i/t 
+A'KRcotnu + 0{R-^) 



, « siuhfl + u)t , 

TTCOtTTI/ / — ; dt 



-R sinh t cosh i/t 



Noting that 



/ ^ , sinh(l + z/)t' 

lim 47rit cot tti/ — tt cot tti/ / at ^-j 

-R-»oo y J-R smh t cosh i/t^ 

/■°° sinh(l - lAt 

= 27rcot7rt/ / at — , 

Jo smh t cosh i/t 
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we have 



1 I r°° , cosht 
Ji = -I tdt 



cotTTZ/ , sinhfl — z/)t 
dt ^ ' 



2 TT^ Jo sinh t cosh^ z/t tt 7o sinh t cosh i^t 
By a similar calculation, starting from 



sinh(l + ^/)(t-7^^) 1 ^ 
c sinh(t — TTz) cosh vt ' 



we obtain J2 = |- Collecting these terms, we obtain (|D.4|) . 



□ 



E Integrals related to the case = 1/2 

Here we supply proofs to the formulas presented in section ^. We retain the 
notation there. 

First let us evaluate the integrals J^i ( p.2|) and J^j ([E.5|) . Set 5^ = e^^ , and 
define 



J. ^7 



1 



(-ir 



nfc=r(-^j + ^k) V^k=^ {Bj - Bk) 

k^3 



{r < j < i < s), 



{r <i<] <s). (E.l) 



Proposition E.l For 1 < k < n and r < i < s we have 



i+l-i 



1/2 



'fei 
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1/2 , 



= 2v^-By'{Y[BA u-Bjr'F,, 



1/2 



j=i 



(E.2) 
(E.3) 

(E.4) 



Proof. Changing the integration variable to A = e"' we have 



-j2r~s+l-i / , \ 1/2 

-i?.'^' n 5. X / ^^^^ 



TT 
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where we have set 

To see ( |E.2| ) it suffices to show that 

This follows from integration of ooki log(— A) along the contour shown in Figure 




Figure 6: The contour for the residue calculus. 



The formula ( [E.4|) is a direct consequence of (|E.2| ). Counting the sum of the 
residues of u^i we find 

j=r j=i 
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This shows the equahty of ( |E.3| ) and ( |E.4| ). □ 
Define 

= 4. + {-ly^'hi, (E.5) 

and denote by U and Jj the column vectors 



Proposition E.2 



1= n 2v^cosh^-^ — X det (J„ J,+i, ■ ■ ■ , J,) . (E.6) 

r<j<k<.s 

Proof. To see this, we write down the known equahty 

i = (i)= E {E^:l---41m.---.Ps)- (E.7) 

■!"' t^s 

For convenience let us introduce the symbols /j(±),/j(±) by setting = 
/j, = /j and /«(—), ) =the empty symbol. Then ( [ti).7[ ) can be written as 

1 = n 2v^cosh^^^ 

r<j<k<.s 

X ^ det (4(-e^),---,J,(-eJ,4(e,),---, J^(e^)) . 

Taking the sum over es,Ss-i, ■ ■ ■ successively and keeping track of the signs, we 
find that the last sum becomes a single determinant det(Jr, ■ ■ ■ , J^). □ 

Proposition E.3 If r < m,l < s then 

{ipmip!) = det {Jr, ■ ■ ■ , Jz-1, /m, Jl+l, ' ' ' , Js) (E.8) 
where D = det(Jr., Jr+i, ■ ■ ■ , Js)- 

Proof. This can be verified in a similar way as in the proof of Proposition ET^ 
As an example, let us take r = 1, s = 4, m = 1, / = 3: 



£2,^4 



Using ( |5.1|) and ( [E.6|) , in the same notation as in Proposition |E.2| , we have 

Di^Pirs) = E^2det(/i( + ),J2(-£2),/3(-),/4(-£4), 
/4(£4),/3(-),/2(£2),/l( + )). 
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The sum over £4 gives J4, Summing further over 62 we obtain 

det(/i, J4, 12, h) - det(/i, J4, h) 
= det(/i, -J2, -''4, A) 
= det(Ji, J2,Ii, J4) 

where we used J2 = h + h and Ji = Ji — Ji . 
In general, consider the case m < I. We have 

D{^rn'ip^) = £m+i- ■ -ei-idetlie), 
where I{e) denotes the matrix consisting of the following array of column vectors: 

//(-), Ii+i{-ei+i), Isi-Ss), Lies), Ii+i{ei+i), 

Summing over Es, Ss-i, ■ ■ ■ we see that the sum combines into a single determinant 

det(t/j-, ■ ■ ■ , Jm—lj Im^ Jm+lj ' ' ' ; ' ' ' ; Js: Im) 

= det( J^, ■ ■ ■ , Jm—ly ( 1) Jm+lj ' ' ' ; Jl — ly -^m; Jl+ly ' ' ' ) <-^s)- 



This shows ( |E.8| ). The other cases are similar. q 



Arguing in a similar manner, one can show in general the following. 

Proposition E.4 Suppose mi < ■ ■ ■ < rrik, li < ■ ■ ■ < 1^, and letr < min(mi, li), 
max(mfc,/fc) < s. Then 

{ipm, ■ ■ ■ Ipm^lpl ■ ■ ■ Ipl) 
= D det (Jj., ■ ■ ■ , J;j_i, Jmi, ■ ■ ■ , </zj.-l, /rrifc; <^is,+l) ■ ■ ■ ) <^s) 

where D = det(Jj., Jr+i, ■ ■ ■ , Js)- In the right hand side Imj is placed at the Ij-th 
slot. 

We omit the details. 

Proof of Proposition \5. 4 Consider the matrix X = {Jr, Jr+i, ■ ■ ■ , Js), and set 
Km = X~^Im- Then Proposition |E.4| states that 



{i^m^ ■ ■ ■ 'ipnik'^l ■■■i^h) = det (ei, ■ ■ ■ , Km^,- ■ ■ , Km^, " " " , e„) 

where the Cj = {Sji)i<i<n denote the unit vectors. It is clear that the right hand 
side is 



rrij )li 



det i(K., 

35 



l<j,i<k 



and that {ipmi'i) = (-^m)/- The proposition follows from this observation. □ 

Proof of Proposition \5. 4 The formula ( |5.5| ) is already known. Let us show 
by taking r = m and s = / in (|E.8|). We wish to compute 



det (t/jnj Jm+l ; ■ ■ ■ ) 1 ? Im 



Substituting 



TT 



Bii' n 



I 



1/2 



1/2 



j=m 



TT 



and using (|E.3| ), we have the following expression for D{iljm4'i)- 

n/2 , 



TT 



Bli' [WbA E (/^. - Pm)G,^ X det X,YZ. 



j—m 



Here X, , Y, Z are the following matrices. 



Y 



/ 1 




1 




1 




—Bm 




-B 


-1 








... (_ 




\n-l 




V 


mm 


Fml-1 


o\ 




































V ■■■ 














diag(2v/^'-sVV.., 27315//^,!). 



1/2 



The matrix Y is upper triangular with diagonal entries F^^ {m < k < I — 1) 
and 1. It is therefore straightforward to compute the determinant. Inserting the 
expressions for Fji and Gji in ( |E.1| ) we obtain the formula ( |5.4| ). 

The case of (ipmi'i) can be shown similarly, using ( [E.4|) . □ 



F Inhomogenoeous Ising model 

In this section, we compute the correlation functions of the inhomogenoeous Ising 
model at the critical temperature. We give an explicit formula for the vacuum 
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expectation values (vac|'?/'i^'?/'„|vac) where ipn^ i'n {n E Z) are the free fermions 
diagonahzing the transfer matrix T{u) of the critical Ising model. The general 



correlation functions are given by the Pfaffians of these matrix elements. In 
the correlation functions for the critical Ising model were given. We have not 
checked the equivalence of our result to theirs except for some simple cases. 

F.l Completely inhomogeneous Hamiltonian 

Consider the transfer matrix of a completely inhomogeneous six-vertex model in 
the infinite volume: 



n+l 











\ 


\ 


\ 





Em 



'^'n'^'n-l / n 



u 



{rn} 



The horizontal line carries the spectral parameter and the vertical lines carry the 
spectral parameters /3„ + u. We assume that /5„ = if \n\ ^ 0. The Boltzmann 

weights RJ^'^l{(3) are given by ( p.3| ) with u 

six- vertex model is equivalent to the critical Ising model (see e.g. [Q). 
Let S be the shift operator 



^. This is the choice in which the 



S 



Then we have 



= ■ ■ ■ Rn+lniPn+1 + u)Rnn-l{Pn + U) 



(F.l) 



where 



/I \ 

m m 
m m 
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and 



The matrix R{I3) can be put in the form 

^(/?) = e^^, X ®a- + ®a+ 

where cr^ = (cr^ ± icr^)/2 and 7 = 7(/3) is related to j3 by 

_^ ^ 1 + isinhf ^ i(l - C^) 

coshf 1 + e l + C ' 

_0 1 — isinh7 

e 2 = = C- 

cosh 7 

As increases from to tt, 7 increases monotonically from to 00. We write 
In — liPn), Cn — cosh 7„ and Sn — sinh 7„. In the below, we assume that 5'„ < 1. 

F.2 Jordan- Wigner transformation 

As usual, we define the Jordan- Wigner transformation 

< = n <' 

fPn^(Jn n 
m<n 

Note that the operators ip* and ipn satisfy the canonical anti-commutation relation 

irm, C]+ = bPm, i^nU = 0, [^^ ^n]+ = ^m.n- (F.2) 

For m, n e Z such that m > n, we set 
Note that 

-^nn—l -^nn— 1) [-^nn— 1) -^m,n] -^mn— 1) [-^nn— 1) -^mn— l] -^mn- (F-3) 

We have 
Proposition F.l 



T(0)-1T(m) = l + ^7i + 0(M') (F.4) 
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where 

Proof. Let us use = to mean an equality modulo u"^. Using ([F.l| ), we have 
T{0)-^T{u) - 1 

n 

X-Rn-ln-2(/5n~l) ' ' ' 

= X/ ' ' ' Rn-ln-2{Pn-l) ^CnHnn-1- 

^ n 

Since Rkk-i{P)~^ = e~^''^^'^-^, the proposition follows from ( |F.3| ). □ 

If we fix /3„'s, the transfer matrices T{u) commute with each other for different 
value of u. Therefore H. also commutes with T{u), and they can be diagonalized 
simultaneously. 

F.3 Diagonalization 

In order to diagonalize the Hamiltonian Ti. we set 

m = E Cue'"' n (1 + Sje-^') n (1 - S,e'')^n, (F.5) 

n j<n-l j>n+l 

Then we have 
Proposition F.2 

[n,m] = -ie'' + e-'')m, (F.6) 

[n,<P*m = {e'' + e~'')<P*ie), (F.7) 

[<P*{e,),<p{d,)U = 11(1 + S,e-''^){1 - S,e^'^)Y.e^'^'^-'-\ (F.8) 

3 k 

Proof. Set z = e*^, and write 

^(^) = E ^"V^"' 

= a^^" n (1+5',^-') n (i-'S',^). 
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Write also, 



A 

We are to prove 



r. 







if m > n; 
if m = n; 



nGZ 



(F.9) 



(F.IO) 



Suppose that /3„ = 0, i.e., C7„ = 1 and Sn = 0, except for M < n < N. If 
m>N + 2 or 'm<M — 2, then (|F.10|) is vahd because 



A 



1 if 72 = m ± 1; 
otherwise, 



and 



M<j<N 



'l + SjZ~^) iin>N + l; 



^"nM<,<^(i-'5,^) ifn<M-i. 



Therefore, ([F.10|) is written in the form 

^(N,M)^{N,M) ^ 

where A^^'*^) = (aI^'^'A and a;(^'^) = (xi^'^A 

matrix A^^''^^ is of the form: 



7V+l>n>Af-l 



j^{N,M) 



—SnCn-1 
SnSn-iCn-2 



Cn 

CnCn-1 
-C]v5jv_iCjv-2 



(F.ll) 
. The 



-^(N-l,M) 



(_l)iV-MS^5jV_,...5„+,CM (-l)^-*^-lCivSjv_i-SAf+lCM 

V (-l)^-''^+i5jv5jv-i-5Af+i5Af {~1)^-''''CnSn-i-Sm+iSm 
( SmSn-i-'-Sm+iCm SmSn-i-'-Sm+iSm \ 

CnSn-i---Sm+iCm CffSf^-i--SM+iSM 

^{N,M+1) 



Cm+iCm 

— z—z~-^ 
Cm 



Cm+2Sm+iSm 
Cm+iSm 
Cm 
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Here 



■j{N-l,M) 



(a{n-i,m) 



N-l>m,n>M- 



^1 — y-"-mn J 



N+l>m,n>M+l 



Similarly, the vector of the form: 



X 



{N,M) 



( {1 + Smz-')---{1 + Snz~')z^+^ \ 
{1 + Smz~^)---{1 + Sn-iz~^)Cmz'' 



V x^N-l,M) X (1 - Snz) 

( x^N,M+l) X (1 + SmZ-') 



J 



\ 



Cm{1-Sm+iz)---{1-Snz)z''' 
V (1 - Smz){1 - Sm+iz) ■ ■ ■ (1 - S^z)z^'-' J 



(F.12) 



Here 



^{N-l,M) 



Af-l>n>A/- 



N+l>n>M+l 



We prove ( [b'.ll| ) by induction on — M. If N = M, we can check directly that 



/ 


-Z-' 


Cn 


Sn 






—z — z~^ 


Cn 


v 




Cn 


—z 



(1 + Snz-')z^+^ 
Cnz^ 



0. 



li N > M, noting that 

— Sn ■ (1 ~l" Snz ^)z^~^^ + Cn ■ Cnz^ ~ (1 — Snz)z^ , 

we can reduce the equality (a^'^^^^x^'^^^A = for N - 1 > n > M - 1 to 
j^^(7v-i,M)^(jv~i,M)j^ _ Q_ Similarly, noting that 

Cm ■ Cmz^ + Sm ■ (1 - Smz)z^'-' = (1 + Smz-^z^ , 
we can reduce the equality =OforA^ + l >n>M+lto 



|^^(Af,M+i)^(iv,M+i)^ _ 0. The proof of (Q is similar. 
Let us prove (F3). We have 



2 A; 
A: ^2^1 



fe-1 



j=—oo j=k+l 
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where Zj = e*^-' (j = 1,2). We assume that Cn = 1 and S'„ = except for 
M < n < N . Then we have 



N+l -N-1^^=m{^ + '5'j^2 ~ SjZi) 



1 - ^2^1 ^ 



+ E 4^r'(i + si) n (1 + ^.^2'')(i - n (1 - s,z,){i + s,^r') 

A:=Af j=M j=k+l 

^ ^Af_l^_M+i n^M(l ~ SjZ2){l + ^) 

^ ^ 1 - 

Write the RHS as 

N 



fceZ j=A// 



Then, F(^'*^)(5m, ■ ■ ■ , 5^; ^i, ^2) belongs to C[5m, ■ ■ ■ , ^iv 

We wish to show F^^) = 0. Let G^^^^''\Sm. ■■■,Sn; ^i, ^2) be the RHS of 
(FU3D . Note that F^^) = G^^'^^) in C{zi, Z2)[SMr ■ ■ , Sn]- Therefore, it is 
enough to show G^^'^^) = as an element of C{zi, ,22) [5*^/, ■ ■ ■ , Sn]. Note that 

G^'^'^^\Sm, ■ ■ ■ , Sn', zi, Z2) = G^^''^^\—Sm,---,—Sn',Z2^,z-^^^) 

= zr''z~,'^'''^G^^^''^\S^,...,SM,Z2.z,). 

It is also easy to show that 

G^'^'^^\Sm, ■ ■ ■ , Sn', zi, Z2) = G^'^'^''^\Scr{M), ■ ■ ■ , S^(^N)', zi, Z2) 

for any permutation a of {M, ■ ■ ■ , A^}. Now G^'^'^'^^Sm, ■ ■ ■ , S^] Zi, Z2) is a poly- 
nomial of degree 2 in Sm- Therefore, in order to show that = 0, it is 
enough to show 

G(^'^n^i,'5Af+i,---,^iv;^i,^2)=0. 
This is shown by induction: we have 

Q^^'^^)(^zi, Sm+1-i ■ ■ ■ , Sn'i Zi, Z2) = {Z2Z1 ^ + Z2Zi)G'^ ' \Sm+i, ■ ■ ■ , •S'at; zi, Z2) 
where 

AT 

G^^'''+'\W,---,5^;^i,^2) = -^2^-'^r^+' {l + S,Z2'){l-S,z,) 

j=M+l 

+ {l-z,z^') E zt'zi'^\l + Sl 



k=M+l 



AT 

X n il + S,z,')il-S,z^) n il-S,Z2)il + S,z^'] 



j=M+l j=k+l 

N 

+ E (i-^,^2)(i + ^,^r^). 

j=M+l 
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Then we can show that 

7^(Af,M+2).„ „ ^ 

Xtj [Dm+2, ■ ■ ■ , ON', Zi, Z2). 

Because G^^'^\Sn; zi, Z2) = 0, we have G^^'^^^^\Sm+i, ■ ■ ■ , S^; Zi, Z2) = by 
induction. 



F.4 Correlation functions 

The vacuum vector |vac) satisfies 

(j){e)\Yac) = if-f<^<f, 0*(e)|vac)=O if f < ^ < f . (F.14) 
Similarly, the dual vacuum (vac| satisfies 

(vac|0(^) = O iff<^<f, (vac|0*(^)=O if -f < ^ < f . (F.15) 

We have also (vac | vac) = 1. Our goal is to compute two point functions (vac|^/'^-?/'n,|vac). 
For this purpose we need 



Proposition F.3 



JO 2,71 



Here An{6) and A* (^) are given by 



a n (1 + S,e-^') n (1 - S,e^') 

j=—oo j=n+l 
^"^^^ " ~ ^ I 1 + Sne-'^ ^ 1 - SnC^' 



n (1 - S,e^') n (1 + Sje-^') 

00 j=n+l 

Proof. From (|F.5|) we have 
dO 

mAnie)- 

Jo C„n?=„+i(l- V) ^l + '^-e-' l-5„e'^ / 27r' 
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Noting that /o^'^e*"^|^ = (5„,o, we can show this is equal to ijjn- The other case is 
similar. 

The two point functions are given as follows. We have obviously 
(vacl^/'mV'nlvac) = (vac I ^/^j^?/^* I vac) = 0. 



Proposition F.4 Suppose that m < n. We have 

{vac\'^*^iljn\y&c) = (-l)'"-"+^(vacKV'm|vac) 

= (-l)™""(vac|V'mdvac) = -{vacl^ni'mW^c) 



Y\m<l<n (^Bj — Bl) 



(F.16) 



Here, we set Bj = e^J . In addition, we have 



(vac I I vac) = (vacjT/'nVnlvac) = ^. 
Proof. Because of (|F.2|), it is enough to compute (vac|-?/'j^^/'„|vac) {m,n G Z) 



Consider the anti-involution 

Ipn ^ i^n, Pn ^ 'Pn, In "Tn, (vac| ^ |vac), 0(6') ^ 4>*{-&). 

Note that the last expression in (|F.16| ) changes the sign by (— 1)"^""^^. Therefore, 
it is enough to prove the equality for (vac|'?/^j^'?/'„|vac). 

First, consider the case m = n. Using ([F.14|) , ([F.15|) and Proposition |F.3| , we 
have 

(vac|'?/'*^/'„,|vac) 

(1 + Sne-''){1-Sne'') ( 1 ,1 ^^M 



2 



d / i + Snc-^' , . \ de 



1 

2 

In general, for m < n, we have 
(vac|V';^V'n|vac) 



\ Cm Vr 1 - 'S'j-e'^ Cn i(n-m)e 

5 1 + ^^e-%.ii,, 1 + 5,e-n + 5„e-^ 27r 



44 



With the change of variable z = — e , the right hand side becomes 

' c„c„n;=^+i(i + 5,z) dz 



- (T\n-m^rnCn_ f ^ + ^ nj=r»+l(l + SjZ) 

Here the branch of log(2; + i)/{z — i) is such that it has the value at 2; = 00. 
The contour C is as in Figure 




Figure 7: The contour C. 



Taking the residues at 2; = Sj {m < j < n) and using (|F.5|) and, in particular, 
the equality 



we have ([F.16|) . 



log : = -Pj - T^h 



J 

j2 nr=m+i(-Bj + Bi) ^ ^ 
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